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Abstract 

We report some new results on the complex-temperature (CT) singularities of g-state Potts 
models on the square lattice. We concentrate on the problematic region Re{a) < (where a = 
e^) in which CT zeros of the partition function are sensitive to finite lattice artifacts. From 
analyses of low-temperature series expansions for 3 < g < 8, we establish the existence, in this 
region, of complex-conjugate CT singularities at which the magnetization and susceptibility 
diverge. From calculations of zeros of the partition function, we obtain evidence consistent 
with the inference that these singularities occur at endpoints Oe, a* of arcs protruding into 
the (complex-temperature extension of the) FM phase. Exponents for these singularities are 
determined; e.g., for g = 3, we find Pe = —0.125(1), consistent with (3^ = —1/8. By duality, 
these results also imply associated arcs extending to the (CT extension of the) symmetric 
PM phase. Analytic expressions are suggested for the positions of some of these singularities; 
e.g., for g = 5, our finding is consistent with the exact value ae,a* = 2(— 1 =i= i). Further 
discussions of complex-temperature phase diagrams are given. 
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1 Introduction and Model 



In this paper, we report some new results on complex-temperature singularities of the g-state 
Potts model |jl|, on the square lattice. The Potts model has been of interest both as an 
example of a particular universality class for critical phenomena and as a model for physical 
phenomena such as the adsorption of certain gases on substrates [0]. However, in contrast 
to the 2D Ising model (equivalent to the q = 2 case) the free energy of the Potts model 
for general q has never been calculated in closed form, even for zero external field(s). Some 
exact results have been established for the model: from a duality relation, the critical point 
separating the disordered, Z^-symmetric high-temperature phase from the low-temperature 
phase with spontaneously broken Zq symmetry and associated nonzero ferromagnetic (FM) 
long-range order is known The free energy and latent heat IQ, and magnetization |Q 
have been calculated exactly by Baxter at this critical point, establishing that the model has 
a continuous, second-order transition for g < 4 and a first-order transition for g > 5. Baxter 
has also shown that although the q = 3 model has no phase with antiferromagnetic (AFM) 
long-range order at any finite temperature, there is an AFM critical point at T = 0. The 
values of the critical exponents (for the range of q where the transition is continuous) have 
been determined 0. Subsequently, further insight into the critical behavior was gained using 
the methods of conformal field theory [0. A review of work up through 1982 was given in 
Ref. i. 

In general, if one knew the exact (zero-field) free energy, one would be able to determine 
the full phase diagram as a function of complex temperature. The idea of generalizing a 
variable on which the free energy depends from real physical values to complex values was 
pioneered by Yang and Lee . These authors considered the generalization of the external 
magnetic field to complex values [§] and proved a celebrated theorem that the complex- 
field zeros of the Ising model partition function lie on the unit circle in the /i plane, where 
H = e~^^'^, pinching the real axis as the temperature T decreases through the critical point 



Tc. Complex-temperature (CT) singularities of Ising models, first considered in Ref. [|10 



were investigated both by means of CT zeros of the partition function |11|-|13| and via their 
effects on low-temperature series expansions [Q. As well as being of historical interest, these 
are relevant here because of the equivalence of the (spin 1/2) Ising model and q = 2 Potts 
model. There is continuing interest in such complexifications because of the deeper insight 
which they give one into the properties of statistical mechanical models (for the Ising model. 



see, e.g., Refs. [0-[^). From general arguments and comparisons with exact solutions for 
2D Ising models with isotropic couplings, one knows that in the thermodynamic limit, CT 
zeros merge together to form curves (including possible line segments) across which the free 
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energy is non-analytic. These curves form the complex-temperature phase boundary (CTPB) 
B of the model. One can define notions of complex-temperature extensions (CTE's) of the 
physical paramagnetic (PM), ferromagnetic (FM) and (if it occurs) antiferromagnetic (AFM) 
phases. In certain cases there are other (labelled "O") complex-temperature phases which 
do not have any overlap with any physical phase. These various CT phases are separated 
by boundaries comprising B. The locus of points making up B may also contain part(s) 
consisting of curves (arcs) or line segments which protrude into and terminate in, certain 
phases. 

There have been several calculations of complex-temperature zeros of the partition func- 



tion for the Potts model on the square lattice |27]-|32]. Since the early calculations for 



q = 3,4, it has been recognized that the zeros show one clear feature: if one uses duality- 
preserving boundary conditions, then in the Re{a) > region (where a = e^; see below 
for notation), these zeros lie on a portion of the unit circle |x| = 1, where x = {a — 1)/^/^ 
p^-|3l|. In passing, we note that in the g ^ oo limit it has been shown (assuming that 
the q ^ oo limit and the thermodynamic limit commute) that the CT zeros lie on the unit 



circle \x\ = 1 |3^, |3^. However, for a given (finite) q, the situation in the Re{a) < region 



has proved to be much more difficult to elucidate. The zeros exhibit considerable scatter 
and, as we shall demonstrate, significant sensitivity to the boundary conditions used for the 
finite lattice calculations, even if one requires these to preserve duality. These facts have 
rendered it problematic to try to make inferences from calculations of zeros on finite lattices 
about the complex-temperature phase boundary B in the thermodynamic limit. In one early 
work it was conjectured that in the thermodynamic limit the zeros lie on the two circles 



I a — 1 1 = y/q and | a + 1 1 = y/A — q for q = 3 and g = 4 (where the second circle degener- 
ates to a point), but shortly thereafter, from a calculation of zeros for the q = 3 model on 
larger lattices, it was concluded that this conjecture was false and the zero distribution 
did not suggest the existence of any simple algebraic expression which would describe this 
distribution. The same conclusion was reached from a calculation of the zeros for the g = 4 
model |9|. 

We have been able to make progress in the problematic Re{a) < region by employing 
a powerful method not hitherto used for this purpose, viz., to combine analyses of low- 
temperature series with calculations of CT zeros of the partition function. We report our 
results here. In our series work, we have taken advantage of the recent calculations of quite 
long low-temperature series for the free energy, magnetization, and susceptibility of the 



square- lattice Potts model for g = 3 up to g = 10 by Briggs, Enting, and Guttmann ^5 



extending earlier calculations (e.g., Ref. for g = 3 and g = 8, see also Ref. [^.) The 
organization of this paper is as follows. In section 2 we define the model and our notation 
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and mention some of the general exact results which are known. In sections 3 and 4 we 
present our new results for the g = 3 and q = 4 Potts model. In section 5 we mention some 
similar results for g > 5 Potts models. Concluding remarks are given in section 6. 

2 Definition of Model and Exact Results 

The (isotropic, nearest-neighbor) g-state Potts model on a lattice A is defined, at a temper- 
ature T, by the partition function 

Z=J: e-^« (1) 

with the Hamiltonian 

^ = -^E^^n<X„, -^E^O<X„ (2) 

{nn') n 

where cr„ = 0, g — 1 are Z^-valued variables on each site n E A; (3 = {ksT)'^; and {nn') 
denotes pairs of nearest-neighbor sites. We use the notation K = pj, h = PH, 

a = z-' = (3) 

and 

and denote the (reduced) free energy per site as / = —^F = limTv^^oo In 2', where Ns 
denotes the number of sites in the lattice. Here we consider the square (sq) lattice. There 
are actually g types of external fields which one may define, favoring the respective values 
(T„ = 0., , ,g — 1; it suffices for our purposes to include only one. The order parameter is 
defined to be 

qM-1 

m = (5) 

g- 1 

where M = (a) = limh^odf /dh. With this definition, m = in the Zg-symmetric, disor- 
dered phase, and m = 1 in the limit of saturated ferromagnetic (FM) long-range order. We 
shall refer to m as the magnetization. Finally, the (reduced, initial) susceptibility is denoted 
as X = l3~^X = dm/dh\h=o- We consider the zero-field model, H = 0, unless otherwise 
stated. 

The general g-state Potts model on the square lattice obeys the duality relation [|l| 

(a-l)(arf-l) =g (6) 
where ad = T^{(i) is the image under the duality map T> of a; 

V{a) = 1 + (7) 
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with = 1 as usuaL In terms of the variable x, this duality relation takes the simple form 

Xd=- (8) 

X 

The critical point at which a phase transition occurs between the high-temperature, sym- 
metric phase and the low-temperature, FM phase, is given by the self-dual point Xc = 1, 
i.e., Qc = I + y/q. One may observe that eq. (||) also has a second self-dual solution at a 
complex-temperature point, 

a = l- ^/q = ai (9) 

Exact expressions for the free and internal energy, latent heat, and magnetization have been 
given by Baxter [|, Q on the critical self-dual curve 

{a-iy = q (10) 

Since the latent heat is zero at Oc for g < 4 [^, the corresponding transition between 
the symmetric and FM phases is continuous. The curve ( ]T0| ) also includes the complex- 
temperature singular point at ae in eq. (^. We note that by the same reasoning, the phase 
boundary associated with the complex-temperature point is also continuous. Exact results 
have also been given by Baxter for the critical manifold |5B 



{a + iy = A-q (11) 

For q = 3, eq. (|lT]) has two solutions: the AFM critical point at a = 0, i.e. T = with 
J < 0, and a complex-temperature point, 

a = -2 = a, (12) 

We observe that these two points are mapped into each other under the duality mapping V: 

V{a = 0) = -2 (13) 

For g = 4, eq. (|Tl|) has only the solution a = —1, implying that for g > 4 the model has no 
AFM long-range order even at T = 0. 

Consider an x Ly planar square lattice G with A^^ = Nq sites (0-cells), A^;, = A^i 
bonds (1-cells) and Np = N2 plaqettes (2-cells). We shall discuss boundary conditions in 
section 3.5 below. The dual of G, which we denote Gd = 'D{G), is defined by associating 
uniquely a (2 — p)-cell of Gd with each p-cell of G, so that (A'o)g = {^2)0^^ (A^i)g = {^1)0^, 
and (A'2)g = (Ao)Gd- Since the planar graph has no handles, Nh = and hence Euler 
characteristic xe = 2(1 — A^^) = 2, it follows from the Euler relation J2'j=oi~^y = Xe 
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that No — Ni + N2 = 2. The duahty relation connecting the partition function Zq on G with 
that on Gd is [H 

ZG(x)=a:^^g^°-i-(i/2)^^ZG,M (14) 

It follows from (|I^) that in the thermodynamic limit, the singularities of that the free energy 
at a point a and its dual image are the same. In particular, 

fsingiq = 3; a = -2) = fsingiq = 3; a = 0) (15) 

In discussing the complex-temperature phase diagram, it is convenient to use the Boltz- 
mann weight z and its inverse a, and the related variable, x. For q = 3, the exact results 
discussed above show that the phase structure for physical temperature (i.e. < a < 00) 
consists of (i) the disordered, Za-symmetric PM phase for < a < 1 + ^/3^, (ii) the FM 
phase for 1 + < a < oo; and (iii) an AFM critical point at a = 0. For g > 4, the 
physical phase structure consists only of the PM phase for < a < 1 + ^/g and the FM 
phase for 1 + y/q < a < 00. One defines the complex-temperature extensions (CTE's) of 
the PM and FM phases by analytically continuing away from the respective segments of the 
positive real a axis. Two rigorous properties are the following. First, because the model 
has a high-temperature series expansion with finite radius of convergence, it follows that the 
CTE of the PM phase occupies a finite neighborhood surrounding the point a = 1. Second, 
it is easy to show that for sufficiently large \a\, one is necessarily in the (CTE of the) FM 
phase. To see this, let a = PaC^^"', then 

K = \na = In Pa + i{ea + 27m) (16) 

where n denotes the Riemann sheet of the logarithm and may be taken equal to zero here. It 
is clear that for sufficiently large \a\ = pa, the angle 6a makes a negligible contribution to K, 
so that (given that d = 2 is above the lower critical dimensionality for the FM transition), 
the system will be in the FM phase. This fact can be seen equivalently as a consequence of 
the fact that the model has a low-temperature expansion with a finite radius of convergence, 
so that there is a finite neighborhood of the origin in the complex z plane where it is in the 
(CTE of the) FM phase. Henceforth, we shall generally refer to the complex-temperature 
extension of the FM phase simply as the FM phase and similarly with the PM phase. 

3 q = 3 Potts Model 

We shall discuss our methods in detail for the g = 3 square-lattice Potts model and then 
proceed to the higher q cases. We begin with analyses of the low-temperature series ex- 
pansions. The series for the partition function, magnetization, and susceptibility have been 
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calculated to order z*^'' |3^. As in our previous studies of complex-temperature singularities 



of various spin models, we have used both dlog Pade and differential approximants to analyze 
the series. We fit the specific heat C, the magnetization m, and the (reduced) susceptibility 
X to the leading singular forms applicable near a generic singular point Zsmg'- 

C r~.il -zlZs^ngY'''-^ (17) 

m ~ (1 - Z/Zsingf'*"^' (18) 

and 

X ~ (1 - Z/ Zsing)"^'''''' (19) 

As usual, the primes indicate that we are approaching this singularity from within the interior 
of the FM phase. We find convincing evidence from our series analyses for singularities at 
two complex-conjugate points, which we denote Zf, and 2;*, at which the magnetization, 
susceptibility, and specific heat are divergent. 



3.1 Magnetization 

In Table 1 we present some diagonal and near-diagonal dlog Pade results for this singularity. 
(We have, of course, also calculated approximants farther from the diagonal.) Our analysis of 
the series for the order parameter indicates complex-conjugate singularities at approximately 
Ze, zl = —0.34 ± 0.29i (the location will be discussed further below) where m diverges with 
the exponent 

Pe = -0.125(1) (20) 
A plausible inference is that the exact value of this exponent is 

Pe = -I (21) 

We comment that although such a divergence in the order parameter is forbidden in usual 
physical phase transitions, it can and does occur at complex-temperature singularities. In- 
deed, in our previous work we have noted several instances where the magnetization diverges 
at CT singularities. For example, exact results show that M diverges (like (1 + 3u)~^^^) at 
the CT point u = —1/3 in the (zero- field) spin 1/2 Ising model on the triangular lattice 
and at u = Ue = —(3 — 2^/2) in the Ising model with (3H = m 12 on the square lattice (like 

(i-v«e)-^/^) m- 

The appearance of the exponent (|20|) at this CT singularity in the square-lattice Potts 
model is intriguing, since —1/8 is not a simple (negative) multiple of any of the physical 
magnetic exponents in the model. This contrasts with the above-mentioned examples from 
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[N/D] 






[13/12] 


-0.340504 + 0.287457i 


-0.1258 


[13/13] 


-0.310120 + 0.287191/ 


-0.1256 


[14/13] 


-0.340382 + 0.2873922 


-0.1250 


[13/14] 


-0.340358 + 0.287473i 


-0.1252 


[14/14] 


-0.340441 + 0.287408i 


-0.1254 


[15/14] 


-0.340428 + 0.287382i 


-0.1252 


[14/15] 


-0.340427 + 0.287365i 


-0.1251 


[15/15] 


-0.340433 + 0.287339i 


-0.1251 


[16/15] 


-0.340436 + 0.287386i 


-0.1253 


[15/16] 


-0.340429 + 0.287368i 


-0.1252 


[16/16] 


-0.340512 + 0.287332i 


-0.1256 


[17/16] 


-0.340405 + 0.287364i 


-0.1250 


[16/17] 


-0.340412 + 0.287339i 


-0.1250 


[17/17] 


-0.340406 + 0.287354i 


-0.1250 


[18/17] 


-0.340416 + 0.287357i 


-0.1251 


[17/18] 


-0.340410 + 0.287350i 


-0.1250 


[18/18] 


-0.340406 + 0.287356i 


-0.1250 


[19/18] 


-0.340407 + 0.287363i 


-0.1250 


[18/19] 


-0.3 10376 + 0.287113/ 


-0.1251 


[19/19] 


-0.340312 + 0.2873252 


-0.1244 


[20/19] 


-0.340260 + 0.287344i 


-0.1241 


[19/20] 


-0.340273 + 0.287344i 


-0.1242 


[20/20] 


-0.340305 + 0.287229i 


-0.1241 


[21/20] 


-0.340285 + 0.287321i 


-0.1242 


[20/21] 


-0.340292 + 0.287319i 


-0.1243 


[21/21] 


-0.340326 + 0.287269i 


-0.1243 


[22/21] 


-0.340711 + 0.286857i 


-0.1246 


[22/22] 


-0.340855 + 0.286878i 


-0.1247 


[23/22] 


-0.340688 + 0.2868882 


-0.1246 


[22/23] 


-0.340982 + 0.287191i 


-0.1260 


[23/23] 


-0.310806 + 0.287311/ 


-0.1262 



Table 1: Values of Ze and (3^ from dlog Pade approximants to low-temperature series for m 
for q — ?,. 
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the 2D Ising model, where, as is clear from the exact solution, the divergent magnetic 
exponents at m = —1/3 on the triangular lattice and aX u = Ue on the square lattice for 
h = iTT/2 are precisely minus the common value of /3 = 1/8 at the physical PM-FM critical 
point. Specifically, for the PM-FM transition in the 2D g = 3 Potts model, the thermal 
and magnetic exponents are yt = 1/v = 1/z/' = 6/5 and yt = 28/15, whence a = a' = 
2 - d/yt = 1/3, S = {d/yh - 1)"^ = 14, /3 = 1/9, 7 = 7' = 13/9, and r] = 4/15 g. One 
recalls that the exponent 1/8 does occur in the set of the conformal weights for the m = 5 
conformal field theory relevant to the 2D q = 3 Potts model, viz., hi 2 = /i4,4 = 1/8, where 
hp^g = [{{m + l)p — mqY — l]/[4m(m + 1)] for p = 1, m — 1, g = 1, and where the 
central charge is given by c = 1 — 6/[m(m + 1)] and has the value c = 4/5 for this case 
0. However, the relation of this to the appearance of the exponent —1/8 at the complex- 
temperature singularities z^,, is obscure, for several reasons. First, as we discussed in Ref. 
|19| , there are violations of basic scaling relations at complex-temperature singularities, so 
that it is not clear how to apply conformal field theory to such singularities (since CFT 
implies, among other things, such scaling relations). Secondly, since the Hamiltonian is 
not real at complex-temperature singularities, it is not obvious why the unitary rational 
conformal series is relevant to such singularities. 



3.2 Susceptibility 

In Table 2 we present our corresponding results from the dlog Pade analysis of the low- 
temperature series for the susceptibility. We have also carried out a similar study with 
(first-order, unbiased) differential approximants, which yields the same value, to within the 
uncertainty. We determine the value of the specific heat exponent at the singularities z^, 
to be 

I'e = 1-14(6) (22) 

where the uncertainty represents a theoretical estimate from the scatter of values among 
different Pade and differential approximants. 

3.3 Specific Heat 

To study the complex-temperature singularities in the specific heat, we have again carried 
out analyses with both Pade and differential approximants. As an illustration, we show in 
Table 3 our results from the latter. Our notation is the same as in our earlier works, e.g. 
Ref. |T^; [L/Mq;Mi\ is the differential approximant to the generic function (j){z) obtained 
as the solution to the ordinary differential equation QQ{z)(f){z) + Qi{z){zd/dz)(j){z) = R{z), 



8 



\NID\ 






[10/ ioj 


n QQ7Q9/1 J_ n 9QnR77'i 

— u.oo ( 0/4 -|- u.zyuo i it 


1 1 9^? 
l.iZO 




n QQQnn8 _l n 08qqqQ/i 
— u.ooyuuo -|- u.zoyoyoz 


1 1 Qn 
i. iyu 


[10/ iuj 


— U.ooyUDZ -|- U.ZooO < oZ 


1 1 7Q 
1.1 / 


[ID/ iDJ 


— u.oooZDy -|- u.zoyoiyz 


1 1 /IQ 




n QQ8981 _L n 98QQ1 81' 
— U.oooZoi + U.Zoyoio? 


1 1 c;n 

l.iOU 


ri s /I Ri 

[i-Oj iOJ 


n QQ81 91 _L n 98Q1 771' 
— u.oooizi + u.zoyi i it 


1 1 Q8 

l.iOO 


fi c; / I 71 
[i5/i7J 


n QQ781 /I _L n 98Q0fi/l^; 


1 19c; 
i.izo 


Fl R /1 71 


n QQ898n _L n 98QQ1 S')' 

— u.oooZou + u.zoyoio^ 


1 1 c;n 

i.iOU 


ri 7 / 1 71 


— Kj.oooA ( U -r U.Zoyoiy? 


1 1 /IQ 

i. i4y 






1 1 1 7 


ri Q / 1 71 

[iy/ i ( J 


n QQ8QQ9 _L n 98Qfi1 Ql' 

— u.ooooyz + u.zoyoio^ 


1 1 fi*? 
1.100 


[1 R / I »1 
[10/ loj 


n '?'?81fi8 _L n 98Q9fi'?^' 


1 1 /I Q 

1.140 


Fl 7 / I 81 


n '?'?899fi J_ n 98QQ1 9^' 

— u.ooozzo -|- u.zoyoizz 


1 1/17 
1.14 ( 


Fl 8 /I 8l 

[io/ ioJ 


n QQ89C^Q _L n 98QQ9C;i' 

— u.ooozoy + u.zoyozo^ 


1 1 /I Q 

1. i4y 


ri Q /I 8l 

[iy/ ioj 


n QQSQQfi _i_ n 9aQQ77.)' 

— U.oooooD -|- U.Zoyo/ 


1.104 


Ton /I si 
[zU/ loj 


— u. 000000 + u.zoyz4zz 


1.101 


Fl 7 /I Ql 
[i ( / iyj 


n QQ8n77 _L n 98Q97Q'i 

— u.oooU( ( + u.zoyz ( o2 


1 1 Q8 


Fl 8 /I Ql 

[io/ iyj 


n QQ8'?'?9 _L n 98QQf;n'i 
— Ki.oooooA -T U.ZoyoOU^ 


1 1 f^Q 
l.lOO 


Fl Q /I Ql 

[iy/ iyj 


n QQ8Q1 1 _L n 98QQQQ') 
— u.ooooii + u.zoyoooz 


1 1 c;9 

1. lOZ 


Fon /I Ql 
[/u/ iyj 


— u. 000000 + u.zoyoyzz 


1.100 


Foi /I ol 
[zi/ lyj 


— u.ooooo/ + u.zoyooyz 


1 1 7n 
1.1 / u 


[18/20] 


-0 338310 + 289320i 


1.151 


[19/20] 


-0.338341 + 0.289364i 


1.154 


[20/20] 


-0.338168 + 0.289333i 


1.145 


[21/20] 


-0.337704 + 0.289629i 


1.126 


[22/20] 


-0.337349 + 0.290843i 


1.111 


[19/21] 


-0.338188 + 0.291209/ 


1.196 


[20/21] 


-0.337205 + 0.289884^ 


1.102 


[21/21] 


-0.340368 + 0.286150i 


1.080 


[20/22] 


-0.336698 + 0.291093i 


1.064 



Table 2: Values of z^. and 7g from dlog Pade approximants to low-temperature series for x 
for g = 3. 
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where Qo, Qi, and R are polynomials of order Mq, Mi, and L. A review of the methods is 
given in Ref. We determine the specific heat exponent to be 



< = 1-0(1) (23) 

As can be seen from Tables 1-3, the magnetization, susceptibility, and specific heat series 
give consistent values for the location of the singularities Ze,z*. Combining these, we infer 
that 

z^, zl = -0.339(2) ± 0.289(2)i (24) 

For our comparison with the plots of zeros of the partition function, it will be convenient 
to re-express this in terms of the a and x variables. We list the results in Table 4 (together 
with positions of the corresponding singularity for higher values of q, to be discussed later). 

In passing, we observe that although we have found violations of scaling relations such 
as a + 2/5 + 7 = 2 and a' + 2/? + 7' = 2 in our previous work at various complex-temperature 
singularities (e.g., Refs. [|1^, |2^), in the present case, we obtain a'^ + 2(3e + 7e = 1-9 ± 0.1 so 
that, to within the uncertainties, this exponent relation is satisfied. 



3.4 Singularities at Dual Images of ae,a* 

A rigorous consequence of the duality of the model is that the free energy also is singular, with 
the same singularity, at the points which are the dual images of the Oe and a*, namely, for the 
central values, P(ae), P(a*) = 0.141±0.462z or equivalent ly, V{xe),V{xl) = -0.496±0.267. 
These points lie in the (CTE)PM phase. Note that since \xc\ = 1 while \xe\ = 0.56, the 
singularities at Xe,x* lie closer to the origin in the x plane than the physical critical point. 



3.5 Connection of Oe, a* Singularities with CT Phase Boundary B 

We would also like to relate these complex-conjugate CT singularities at Ze,z*, or equiv- 
alently, in the complex a plane, at ae,a*, to the complex-temperature phase boundary B. 
From our previous studies on CT singularities |2^, |2^, |2^, we formulated a conjecture that 
whenever an arc or line segment of the phase boundary protrudes into, and ends in, the FM 
phase, there is a divergence in M at the endpoint of this arc. We proved that this divergence 
in M implies also a divergence in x the same endpoint [^]. Besides the exact results 
alluded to above which exhibit this behavior, our calculations of partition function zeros for 
the 2D higher-spin square-lattice Ising model P3|, in conjunction with the series analyses of 
Jensen, Guttmann, and Enting are consistent with the conjecture. A natural extension 
of the conjecture is that the divergences which we have found in m at these CT points in the 
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[L/Mo] Mi\ 




"e 


[8/16; 14] 


-0.3407788 + 0.2885827z 


1.072 


[8/16; 15] 


-0.3407563 + 0.2878427z 


1.052 


[8/16; 16] 


-0.3395542 + 0.2877949i 


0.997 


[8/16; 17] 


-0.3399969 + 0.2874121? 


0.986 


[8/17; 15] 


-0.3393125 + 0.2891874i 


0.960 


[8/17; 16] 


-0.3391475 + 0.2886514? 


0.967 


[10/15 


13] 


-0.3105761 + 0.2887910/ 


1.019 


[10/15 


14] 


-0.3405514 + 0.2889317? 


1.044 


[10/15 


15] 


-0.3400580 + 0.2896685? 


0.987 


[10/15 


16] 


-0.3392558 + 0.2894922? 


0.935 


[10/16 


14] 


-0.3396813 + 0.2885613? 


0.989 


[10/16 


15] 


-0.3397607 + 0.2885435? 


0.994 


[12/14 


12] 


-0.3401737 + 0.2893914? 


1.001 


[12/14 


13] 


-0.3405823 + 0.2891851? 


1.049 


[12/14 


14] 


-0.3406363 + 0.2897800? 


1.059 


[12/14 


15] 


-0.3395286 + 0.2896329? 


0.934 


[12/15 


13] 


-0.3407441 + 0.2896263? 


1.071 


[12/15 


14] 


-0.3404562 + 0.2891536? 


1.039 


[14/13 


11] 


-0.3399260 + 0.2888016? 


0.989 


[14/13 


12] 


-0.3399969 + 0.2885408? 


1.001 


[14/13 


13] 


-0.3395575 + 0.2895497? 


0.919 


[14/13 


14] 


-0.3395242 + 0.2893949? 


0.915 


[14/14 


12] 


-0.3405031 + 0.2889864? 


1.030 


[14/14 


13] 


-0.3395982 + 0.2892116? 


0.942 


[16/12 


10] 


-0.3394602 + 0.2889824? 


0.955 


[16/12 


11] 


-0.3391883 + 0.2891054? 


0.929 


[16/12 


12] 


-0.3394838 + 0.2892739? 


0.939 


[16/12 


13] 


-0.3393726 + 0.2892381? 


0.934 


[16/13 


11] 


-0.3394488 + 0.2892846? 


0.937 


[16/13; 12] 


— U.Ooy4oD4 + U.ZoyoU40? 


u.yoo 


[18/11; 9] 


-0.3404461 + 0.2882994i 


1.032 


[18/11; 10] 


-0.3398013 + 0.2888885i 


0.979 


[18/11; 11] 


-0.3395252 + 0.2888272i 


0.962 


[18/11; 12] 


-0.3396986 + 0.2891254? 


0.955 


[18/12; 10] 


-0.3394656 + 0.2892506? 


0.939 


[18/12; 11] 


-0.3394441 + 0.2891136? 


0.945 



Table 3: Values of Zg and a'^ from 
(reduced) specific heat, C /{ksK'^), 



differential approximants 
for g = 3.11 



to low-temperature series for 



Q 








3 


-0.339(2) ±0.289(2)i 


-1.71(1) T 1.46(l)i 


-1.56(1) T0.841(3)i 


4 


-0.288(2) ±0.270(2)i 


-1.85(1) T 1.73(1) 


-1.42(1) T0.866(4)i 


5 


-0.251(2) ±0.251(2)i 


-1.99(1) T 1.99(l)i 


-1.34(1) T0.891(4)i 



Table 4: Values of Zg and, correspondingly, and Xe from analyses of low-temperature series 
for magnetization, susceptibility, and specific heat. See text for further details. 



square-lattice g-state Potts model indicate that these points are endpoints of arcs of points 
where / is non-analytic, i.e. arcs on the CT phase boundary B. 

In order to test this conjecture here, we have carried out new calculations of complex- 
temperature zeros of the partition function for several different finite lattices, using transfer 
matrix methods, as in our earlier study of partition function zeros for the 2D higher-spin Ising 
model At appropriate points we shall make comparison with previous computations of 
zeros for the square-lattice Potts model [^-[^. It is desirable to, and we shall, restrict 
ourselves to lattices with duality-preserving boundary conditions (DBC's) 0. These 

guarantee that T>{G) = G i.e. the dual of a (finite) lattice G is (isomorphic to) the original 
lattice. To discuss these, we recall eq. (p^ ) and the associated definitions. Note that 2D 
periodic, i.e. toroidal, boundary conditions do not preserve duality for a finite lattice. The 
boundary conditions are not uniquely specified by the condition that they preserve duality. 
One type, which we label as type 1, was discussed in Ref. ||HT[] (see their Fig. 1). We shall 
need a straightforward generalization of it to the case of an x Ly lattice with ^ Ly, 
and we describe this as follows. Let the lattice be oriented with the x and y directions being 
horizontal and vertical, respectively. Let all of the sites on the upper and right-hand edges, 
including the corners, connect along directions outward from the lattice to a common site 
adjoined to this lattice (so that the upper right corner connects to this adjoined point via 
bonds in both the x and the y directions), while all of the sites on the lower and left-hand 
edges, excluding the previously mentioned corners, have free boundary conditions. For the 
dual lattice, the special adjoined point may be taken to lie to the lower left of the lattice. 
In Ref. ||31| it was noted that from (unpublished) calculations of zeros with other DBC's, 
they obtained agreement with their conclusions from type 1 DBC's that for Re{x) > the 



zeros lie on the circle \x\ = 1 |31|. We have also used a set of DBC's different from type 1, 



which we denote as type 2 |34]. For these, let the x Ly lattice have periodic boundary 
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conditions (PBC's) in the x direction, so that the lattice may be pictured as a cyhnder with 
its axis oriented vertically. Now connect all of the sites on the upper edge of the cylinder 
to a special point adjoined to the lattice, and let all of the sites on the lower edge of the 
cylinder have free boundary conditions. The dual lattice is constructed in the usual way, 
assigning sites to each 2-cell of the original graph, and adjoining the special point below the 
cylinder with the stipulation that the points on the upper edge of the dual lattice have free 
boundary conditions and the points on the lower edge connect to the adjoined point. For 
both type 1 and type 2 DBC's, Nq = L^Ly + 1 = N2 and A^i = 2LxLy so that the prefactor 
^Afo-i-(i/2)iVi = 1 in (|||)_ Note also that both type 1 and type 2 DBC's force some 
sites to have coordination number 3 rather than 4, and the adjoined point has coordination 
number + Ly for type 1 DBC's and for type 2 DBC's. This is in contrast to periodic 
boundary conditions, which violate duality but maintain equal coordination number for all 
lattice sites. 



A third type of DBC has recently been suggested to us by F. Y. Wu |24|; it can be 
defined as follows, and will be denoted as DBC type 3: consider an x Ly lattice with 
Lx = L + 1, Ly = L. Let all of the sites on the longer (horizontal) upper edge of the lattice 
be connected to a special adjoined point, via L + 1 bonds, and similarly, let all of the sites on 
the lower edge of the lattice be connected by L + 1 bonds to a second adjoined point. Finally, 
connect the two adjoined points by a single bond and let the sites on the vertical edges of 
the lattice have free boundary conditions (in the outward horizontal directions). This lattice 
has No = N2 = L{L + 1) + 2 and A^i = 2{L'^ + L + 1) = 2{Nq - 1) (so again, qNo-i-(i/2)N^ ^ ^ 
in ([T^)). Evidently, type 3 DBC's share greater similarity with type 1 than type 2 since in 
types 1 and 3 no subsets of edge sites have periodic boundary conditions, while in type 2 
the sites on the vertical edges do have PBC's. For a given value of q, the patterns of zeros 
which we have obtained with type 3 DBC's are indeed, similar to those with type 1. 

In Figs. 1(a) and 1(b) we present calculations of complex-temperature zeros, in the 
complex a plane, of the partition function for the g = 3 Potts model on 8 x 10 lattices with 
duality preserving boundary conditions of types 1 and 2, respectively. The positions of the 
singular points marked with small circles on both of these plots. A comparison of 

these plots gives a quantitative measure of how the positions of the zeros can vary for different 
boundary conditions, and specifically for those which maintain duality. This extends previous 
published plots [BD, RT|, which showed that the pattern of zeros differs significantly when 



one uses duality-preserving, as opposed to duality-violating, boundary conditions. These 
comparisons demonstrated that once one specializes to duality-preserving BC's, the zeros in 
the Re{a) > region lie nicely on the circle |a — 1| = ^/q, whereas they lie close to, but not 
exactly on, this circle for duality-violating BC's. Among previously published calculations 
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Figure 1: Plot of zeros of Z. in the a plane, for the square-lattice g = 3 Potts model on an 
8 X 10 lattice with duality-preserving boundary conditions (DBC's) of (a) type 1; (b) type 
2; and (c) on a 7 x 6 lattice with type 3 DBC's. The singularities marked with o. 
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with duality-preserving BC's, there is one, given as Fig. 11.1 in Ref. |3^, on a lattice as large 
as the one that we use, and the pattern of zeros found there is very close to the pattern in our 
Fig. 1(b). In our figure, one can discern two outer-lying complex-conjugate arcs of zeros in 
the "northwest" and "southwest" quadrants, and the points Oe, a* lie near the ends of these 
arcs. In Fig. 1(a) there is more scatter among the zeros, but nevertheless, the points ae-,a* 
lie at the ends of subsets of zeros which can be associated with arcs. We have also performed 
analogous calculations of zeros on smaller lattices with similar results. For comparison, in 
Fig. 1(c) we show an exploratory calculation with type 3 DBC's. Taking into account that 
the lattice for Fig. 1(c) is somewhat smaller than that for Figs. l(a,b), one sees that the 
distribution of zeros is similar to that in Fig. 1(a) with type 1 DBC's, especially in the region 
near ag, a* and their dual images. From all of these calculations, we thus conclude that the 
observed patterns of zeros are consistent with the hypothesis that in the thermodynamic 
limit, the points the ends of arcs contained in the complex-temperature boundary 

B which protrude into the FM phase. By duality, this is equivalent to the statement that 
the singularities at V{ae), ^(^^e) endpoints of arcs contined in B which protrude into 

the PM phase. 

3.6 Singularities at a = and a = —2 

The exact solution (^) by Baxter shows that the point a = 0, i.e. T = for J < 0, is the 



antiferromagnetic critical point. By the duality property (|1^), it follows that the free energy 
is also singular, with the same singularity, at the dual to this point, namely, a = —2 = 
given in eq. ([l^). From renormalization group methods, a mapping to a critical 6- vertex 
model, and studies of correlation functions, it has been concluded that at the AFM critical 
point is an essential singularity, with an essential zero in the free energy and an exponential 
divergence in the correlation length as i^' ^ — oo, i.e., a \ [^]. Thus, if one assigns 
algebraic exponents a and u for this AFM critical point, then a = — oo, z/ = oo. By duality, 
the same singularity in the free energy occurs at the dual point a = —2. 

We have addressed two questions concerning these singularities at a = and —2: (i) 
how do they connect with the complex-temperature phase boundary B, (ii) how well do the 
low-temperature series detect the CT singular point at a = —2? For question (i), we first 
recall that the density g of CT zeros along the curves comprising B in the vicinity of a generic 



singular point behaves as W^ 



(7~ |a-a,|^-"^ (25) 

where denotes a singular point and as (a^) denotes the corresponding specific heat expo- 
nent for the approach to from within the CTE PM (FM) phase. Since a = — oo at a = 
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and hence, by duality, a'^ = — oo at = —2, there is a strong reduction in the density of CT 
zeros as one approaches these respective points a = and a = —2 along B. This is consistent 
with what is observed with CT zeros calculated on finite lattices; one sees clear arcs of zeros 
and, e.g. for type 1 duality-preserving boundary conditions, these track toward the respec- 
tive points a = and a = —2, ending some distance away from these points. (For type 2 
DBC's, one also observes a slight curling tendency among the last few points on the curves.) 
From this tracking of the zeros toward a = (and hence, by duality, toward a = —2), we 
infer that the AFM critical point a = lies on a portion of B which connects with the 
continuation of the curves lying on (either part or all of) the unit circle |a — 1| = in such 
a way as to bound completely the complex-temperature extension of the PM phase. This is 
reminiscent of another model which has no AFM long-range order at any finite temperature, 
but an AFM critical point at T = namely, the (isotropic, spin 1/2) Ising model on the 
triangular lattice. In that case, in terms of the variable u = e~^^, the complex-temperature 
phase boundary B consists of the union of the circle \u + 1/3| = 2/3 and the semi-infinite 
line segment running along the negative real u axis from —1/3 to — oo, or, equivalently, in 
the variable analogous to a, it consists of the union of the circle \u^^ — 1| = 2 and the 
line segment —3 < u^^ < 0. The AFM critical point at u^^ = forms the right-hand end 
of this line segment and is connected to the rest of S by it. It is interesting to contrast this 
with situation in a model which is disordered and noncritical at = — oo, i.e. at = 0. 
An example is provided by the (isotropic, spin 1/2) Ising model on the 3 ■ 6 ■ 3 • 6 (Kagome) 
lattice |E3|, WM. Since the point = is noncritical, it must be true that this point is not 



connected to B, and, indeed, one finds [18, 21 1 that is lies in the interior of the symmetric 
PM phase. The same behavior is found for the 3 ■ 12^ lattice ||2l[| . 

Returning to the q = 3 Potts model, the property that the AFM critical point at a = 
is connected to the rest of B as described above imphes, by duality, that the singular point 
at a = —2 lies on the dual image of the above-mentioned curve, which therefore encloses a 
complex-temperature phase and separates it completely from the (CTE of the) FM phase. 
Given the scatter of the zeros, it is not possible to make a very reliable inference for where 
these curves intersect the continuation of the circle |a — 1| = \/3 in the "northwest" and 
"southwest" quadrants of the complex a plane. We do remark that the zeros are consistent 
with the possibility that these intersection points are a = e='=^'^*/^. 

As regards the analysis of the low-temperature series, we find that these series are able to 
locate the singular point at = —2, i.e., Zs = —1/2, but not very accurately. Representative 
results for diagonal dlog Pade approximants are given in Table 5, again based on the series 
calculated to 0{z'^'^) in Ref. [^. We do not show the results of the dlog Pade approximants 
[N / D] with N = D < 13 because these did not locate the z = —1/2 singularity with accept- 
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\N/D\ 


Zs 


„ / 

7. 


[14/14] 


-0.478716 


1.66 


[ib/ibj 


— u.4yuooo 


Z.Ol 


[16/16] 


-0.487890 


2.12 


[17/17] 


-0.487889 


2.12 


[18/18] 


-0.487906 


2.12 


[19/19] 


-0.487844 


2.11 


[20/20] 


-0.487995 


2.12 


[21/21] 


-0.485409 


1.99 



Table 5: Values of Zg and 7^ from diagonal dlog Pade approximants to low-temperature series 
for X for g = 3. 

able accuracy (e.g., the [12/12] and [13/13] approximants gave —0.607830 and —0.572062, 
respectively). As is evident from Table 5, the series yield values of the susceptibility exponent 
7^ ~ 2. Our previous work [|l^ has shown that there are subtleties in applying exponent 
relations like u'd = 2 — a' and 7' = i/'(2 — rj) at complex-temperature singularities. However, 
it is of interest to note that, since a'^ = —00, if the hyperscaling relation u^d = 2 — a'^ is valid, 
then u'g = 00 and hence, assuming that the relation 7^ = i^^(2 — rjs) holds at Zg, it would 
follow that 7^ = 00, i.e., x would have an exponential divergence at this point. Of course, the 
series analysis cannot yield 7^ = 00, but it does produce a rather large value. For reference, 
one may recall that in the case of the 2D 0(2) model, while the Kosterlitz-Thouless theory 
implies that x diverges exponentially as T \ T^, so that 7 = 00, the earlier series analyses 
gave a value of roughly 7 = 3 [^] . 

3.7 Singularity at a = 1 - S^/^ 

The Baxter solution (|1^) shows that the free energy is also singular at the self-dual complex- 
temperature point Ui given in eq. (||) for q = 3, viz., a = 1 — v^. There is considerable scatter 
of zeros in the vicinity of this point (less for Fig. 1(b) or Fig. 11.1 of Ref. [^ than in Fig. 



1(b)); however, the zeros are consistent with the inference that this singular point lies on a 
segment of the circle [a — 1| = V^- As one approaches the regions near the intersection points 
discussed above, the scatter of zeros becomes too great to draw a firm conclusion about this 
part of B. Nevertheless, we are able to infer that the point ae is completely separated from 
the FM phase by portions of the CT phase boundary B. To see this, assume the contrary. 



17 



i.e. that one can analytically continue from z = to zi = l/a^. First, this would contradict 
the property that the singular point at a = —2 lies on a portion of B connecting it to the rest 
of B (the dual image of the curve connecting the physical AFM critical point to the rest of 
B). Second, if were not completely separated from the FM phase, then one should be able 
to detect this singular point with the very long low-temperature series available. However, 
we found no evidence for a singularity at this point from our analysis of these series. The 
obvious conclusion is that ai lies in a region beyond the applicability of these series, i.e., 
beyond the border of the (CTE)FM phase. 



3.8 A Comment on dim{B) for Re{a) < 

We comment here on another feature of the pattern of zeros. For the 2D spin 1/2 Ising model 
with isotropic spin-spin couplings, one knows from exact solutions that on most lattices, the 
zeros merge to form a 1-dimensional algebraic variety, i.e. the CT phase boundary B. Even 
for isotropic couplings, there is a heteropolygonal lattice, i.e., the 4 ■ 8^ lattice, for which 
this is not the case; the locus of points where the free energy is non-analytic forms a 2- 
dimensional algebraic variety [^. Moreover, for non-isotropic spin couplings, this is also 
true, even on the square lattice [^. In both cases it is easy to see why this is true (see 
section 6 of Ref. |^). The zeros of the 2D Ising model for higher spin values also appear to 
approach curves as the lattice size gets large P3 . For the 2D g-state Potts model, the zeros 



in the Re{a) > half-plane lie on a 1-dimensional curve, i.e., part of a circle p8[-||33|. In the 



Re{a) < region, we are not aware of any proof of this. However, we can observe that in 
the known cases with exact solutions, in the thermodynamic limit, the zeros either (i) form 
curves, or (ii) areas, but not both curves and areas. Thus, given that the zeros for Re{a) > 
do form a curve, one would have a qualitatively new situation not previously encountered if 
some of the zeros in the Re{a) < did form areas. 



4 g = 4 Potts Model 

We have carried out exploratory analyses for higher-g Potts models on the square lattice and 
have found evidence for singularities analogous to Ze, z* in each of the cases studied. We 
begin with the q = 4 model. Like the g = 2, 3 cases, this model has a continuous, second- 
order PM-FM transition. However, from the exact solution by Baxter on the manifold (|TT|), 
it follows that the g = 4 model does not have any AFM critical point. Note that the 
complex-temperature point a = — 1 lies on the manifold (|n|). 

We have analyzed the low-temperature series for the magnetization, susceptibility, and 
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specific heat as before. Tliese yield a consistent indication of singularities at the points 
listed in Table 4. (The values of Zf. from the x series exhibit somewhat greater scatter 
than those from the m and C series.) It is interesting that, to within the uncertainty, our 
numerical value for Im{al) can be fit by the exact expression Im{al) = for g = 4. 
This and a similarly intriguing result which we find for q = 5 suggest that there may be 
simple algebraic formulas for the singularities ctg , • view of the logarithmic confluent 
singularities which are known to occur at the physical PM-FM transition in this model 



47| , it would be useful in future work to carry out a more sophisticated analysis of the 



series including such confluent singularities also at the complex-temperature singularities; 
however, in the present exploratory work we have not done this. From our study of the 
low-temperature series for m, we find the exponent Pe = —0.12(1), so that, as before, the 
magnetization diverges at Ze,z*. This rigorously implies that x also diverges at these 
points. 

In order to see how these singularities connect to the complex-temperature phase bound- 
ary B for the model, we have carried out a calculation of CT zeros of the partition function 
on various lattices with duality-preserving boundary conditions. In Figs. 2(a) and 2(b) we 
show calculations of these zeros, in the complex a plane, for a 6 x 6 lattice with type 1 and 
type 2 DBC's, respectively. Previously published calculations of zeros for this model include 
plots for strips (L^ x 32 for = 4, 6, and 8; and 10 x 16) and a plot for a 4 x 4 lattice 
with type 2 DBC's. We use symmetric lattices since in taking the thermodynamic limit on an 
Lx X Ly lattice, if L^/Ly deviates strongly from unity, the results can involve 1-dimensional 
artifacts. For type 1 DBC's (Fig. 2(a)) one sees a clear indication of a complex-conjugate 
pair of arcs of zeros in the "northwest" and "southwest" quadrants, with the singularities 
at Oe, a* forming the endpoints of these arcs. These arcs are not as clear with type 2 DBC's 
((Fig. 2(b)), but again, the points ae,a* lie at the ends of subsets of zeros which can be 
associated with arcs. A calculation with type 3 DBC's is included as Fig. 2(c). The re- 
sultant pattern is very similar to that in Fig. 2(a) with type 1 DBC's, to an even greater 
extent than in the g = 3 model. The patterns of zeros in Figs. 2(a)-2(c) are all in good 
agreement with the above-mentioned property that the model is not critical even at a = 0, 
i.e., that it has no AFM critical point even at T = 0. For the plot in Fig. 2(b) with type 
2 DBC's, the zeros are consistent with the expectation from the Baxter solution ( |Tll) that 
the complex-temperature phase boundary B passes through the point a = —1 (equivalently 
for this q = 4 case, x = —1). In contrast, for the plots in Figs. 2(a) and 2(c), with type 
1 and type 3 DBC's, respectively, there are no zeros near to, or easily extrapolated toward, 
a = — 1. A clarification of the situation in the vicinity of this point merits further study on 
larger lattices. 
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Figure 2: Plot of zeros of Z, in the a plane, for the square-lattice g = 4 Potts model on a 
6x6 lattice with DBC of (a) type 1; (b) type 2; and (c) on a 6 x 5 lattice with DBC's of 
type 3. Notation is as in Fig. 1. 
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5 Higher-Q' Potts Model 



We proceed to square-lattice Potts models with g > 5, for which the physical PM-FM phase 
transition is first-order. Of course, a series analysis does not, in general, yield an accu- 
rate determination of the location of the phase transition point for a first-order transition. 
However, there is no obvious reason why this should be a drawback for our study of the 
complex-temperature arc endpoint singularities Ze-,zl, since these yield a strong signal in the 
form of a divergent magnetization. In Table 6 we list the results for z^. and corresponding 
exponent (3e from diagonal dlog Pade approximants to the low-temperature series for m. 
We obtain similar results for z^ from the susceptibility and specific heat series. From these 
diagonal (and near-diagonal) dlog Pade approximants, we we obtain the value for Ze listed 
in Table 4. This is an intriguing result, since, to within the uncertainty, our determination 
is consistent with the following exact analytic formula: 

ae,a* = 2(-l ^z) for g = 5 (26) 

Of course, in the absence of an exact solution of the model, we cannot exclude the possi- 
bility that this agreement is fortuitous, but it motivates one to think further about simple 
analytic expressions for the location of the arc endpoints which we have discovered. From 
our series analysis, we find that the magnetization diverges at these points with exponent 



/3e = —0.11(1). As noted above, this implies ||2^ that x must also diverge at these points 



and our analysis of the low-temperature series for x yields the exponent 7^ = 1.2(1). 

We find the same generic features for all of the q values that we have analyzed, viz., 
complex-conjugate singularities at points ^e, zl. These thus appear to be a general feature of 
the g-state Potts model on the square lattice for g > 3. Besides the case g = 5, we have made 
exploratory studies of the low-temperature expansions for the cases g = 6, 7, and 8. We find 
z^, zl = -0.23 ± 0.24i, -0.21 ± 23z, and z^ = -0.20 ± 22i for g = 6, 7, and 8, respectively. 
Thus, the magnitude \ze\ (|ae|) decreases (increases) as g increases. Note that if, in the x 
plane, the complex-conjugate arcs retract toward their respective points of origin and finally 
disappear in the q ^ 00 limit, as is required in order for the complex-temperature phase 
boundary B to reduce to the unit circle |x| = 1 in this limit, it is necessary that for large g, 
\Re{ae)\ and \Im{ae)\ grow like ^g. We find that Pe increases from ~ —0.11 for g = 6 to 
~ -0.10 for g = 8. 

In Figs. 3 and 4 we show our calculation of CT zeros of Z for g = 5 and g = 6 on a 
6x6 lattice with type 2 DBC's. (We have checked that other DBC's give results which also 
support our conclusions.) Although there is considerable scatter of zeros in the Re{a) < 
region, the pattern is again consistent with the conclusion that the singularities at a^, a* lie 
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[16/16] 


-0 251668 + 251140i 


-0.1199 


[17/16] 


-0.251631 + 0.251079Z 


-0.1191 


[16/17] 


-0.251667 + 0.251138i 


-0.1199 


[17/17] 


-0.251710 + 0.251167i 


-0.1205 


[18/17] 


-0.251396 + 0.251142i 


-0.1181 


[17/18] 


-0.251644 + 0.251144i 


-0.1198 


[18/18] 


-0.252009 + 0.251663i 


-0.1282 


[19/18] 


-0.251915 + 0.251720i 


-0.1267 


[18/19] 


-0.251910 + 0.251717i 


-0.1266 


[19/19] 


-0.251929 + 0.251724i 


-0.1270 



Table 6: Values of Ze and (3e from dlog Pade approximants to low-temperature series for m, 
for g = 5. 
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Figure 3: Plot of zeros of Z, in the a plane, for the square-lattice q — 5 Potts model on a 
6x6 lattice with DBC of type 2. Notation is as in Fig. 1. 



-2.5-2.0-1.5-1.0-0.5 0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 
Re(a) 



Figure 4: Plot of zeros of Z, in the a plane, for the square-lattice q — 6 Potts model on a 
6x6 lattice with DBC of type 2. Notation is as in Fig. 1. 
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at the end of arcs of zeros, so that in the thermodynamic hmit these points are endpoints of 
arcs of singularities which connect with the rest of the complex-temperature phase boundary 
B. 



6 Conclusions 

In summary, we have used analyses of low-temperature series expansions to study complex- 
temperature singularities in the square-lattice g-state Potts model. We have found singu- 
larities at complex-conjugate pairs of points and, by means of comparison with patterns of 
partition function zeros, have obtained support for the inference that in the thermodynamic 
limit these are endpoints of arcs lying on the complex-temperature phase boundary B. At 
these points, the magnetization diverges, in agreement with an earlier conjecture which we 
had formulated. This guarantees that the susceptibility also diverges at these points, and 
our series analyses are in accord with this. Our work includes several intriguing findings, 
including the likely exponent value /3e = —1/8 for g = 3 and an inference of an exact formula 
(p6|) for the endpoint singularities for g = 5. From the duality of the model, it follows that 
these arcs protruding into the complex-temperature extension of the FM phase are accom- 
panied by their dual images, i.e. arcs protruding into the CT extension of the PM phase. 
Our new results further elucidate the complex-temperature phase diagrams of square-lattice 
Potts models. 



We thank Prof. F. Y. Wu for the helpful private communication This research was 
supported in part by the NSF grant PHY-93-09888. 
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